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Abstract 

We propose a phenomenological two component model describing the decay amplitude of the process 
—t Stt, whose rate has been found surprisingly large. The first component is a constant 
background Fnr, and the second one is a Breit-Wigner type amplitude associated to a quasi 
two body /o(980) 7r+ state. We show that it is possible to reproduce the observed rate for I?+ 
TT+TT+TT^ as well as the two other measured branching ratios for the non resonant part and the 
resonant /o7r+ one, with a common parameter -Fjvi?- 

Predictions are given for the Df — > tt'^tt'^tt^ rates, as well as for the various and tt^ , or tt'^ and 
TT"*" energy distributions for the two possible final states. 
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I Introduction 



We are interested, in this paper, in the decay of the meson into three pions. We observe that 
both flavoured constituents c and s of the meson are absent in the final state. As a consequence, 
the decay mode ^ 3 vr can be described neither by a spectator, nor by a colour suppressed, nor 
by a penguin diagram, but only by a annihilation where the virtual time-like decays into three 
pions. 

It is generally expected that W annihilation amplitudes are small compared to spectator and 
colour suppressed ones. Surprisingly experimental data do not seem to support this belief. As an 
illustration, let us consider the three following decay modes of the D'^ meson involving the same 
Cabibbo, Kobayashi, Maskawa (CKM) factor V*gVud and compare their branching ratios [Q] : 

i) Spectator ; 

Br{Dt ^ (/) + 7r+) = (3.5 ± 0.4) % , 

ii) Colour suppressed ; 

Br{D+ -^K° + K+) = (3.5 ± 0.7) % , 

iii) W annihilation ; 

Sr(L>+ ^ 7r+ + 7r+ + 7r") = (1.35 ±0.31) % . 

It is therefore interesting to understand the origin for such a large value of the — > Stt branching 
ratio. 

The transition of the virtual W into three pions occurs also in the r decay mode r+ — > F,- + 37r. 
However, in this case, the full axial vector weak current contributes, and the amplitude depends on 
three independent structure functions. For the decay Df Sir, only one structure function associated 
to the divergence of the axial vector current is involved, because here the W is in the = 0~ state, 
like the Dg meson. Therefore the study of Dg 3tt decay is the cleanest way to isolate the structure 
function F4 of references ^. 

In Section II, we develop the formalism for D'^ 3tt decay. We study the kinematics and the 
three pion phase space. Experimental data as quoted in Reference Q are listed. 

By inspection of these experimental data, we observe that the quasi two body state p°tt~^ assumed 
in ^ to be responsible of the structure function F4 is negligible compared to the quasi two body state 
/o(980)7r+. Moreover, since the non resonant rate is large, we propose a model in Section III where 
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the decay amplitude for Df Sir has two components : the first one is associated to a background 
non resonant three pions, and the second one corresponds to a quasi two body /o tt"*" final state. 

We compare our model with experiment in Section IV for the final state tt+tt+tt" where data on 
rates are available and we make predictions for the tt'^ and 7r~ energy distributions. 

For the final state tt^tt^tt'^, using isospin analysis, we make predictions in Section V for the rates 
and the 7r° and tt'^ energy distributions. 

Our model is consistent with the rate measurements. However the serious test would be the 
observation of the various tt meson energy distributions and ultimately, if statistics is copious enough, 
of the Dalitz plots. 



II Generalities and Kinematics 

We study the decay of a meson of energy momentum P into three pions of energy momenta 
Pi,P2,P3 with the relation P = pi +P2 +Pz- We introduce the Mandelstam variables si, S2, ^3 and the 
TT meson energies Ei,E2, in the rest frame. Neglecting the mass difference between charged 
and neutral pions, we get 

si = {P2 + Psf = {P - Pif = mjj^ +ml- 2mD,Ei , 

52 = {p3+Pif = {P-P2f = ml^+ml-2mD,E2 , (1) 

53 = (Pi + P2f = (P - Psf = +ml- 2mD,E3 . 

Energy momentum conservation implies the relations 

si + S2 + S3 = m|)^ + 3ml > Ei + E2 + E3 = niD, ■ (2) 

The double differential distribution is given by 

dr = —^ |<37r|T|L>+ >|^ dEx dE2 , (3) 

64 TT"* 

and the transition matrix element < ?>tt\T\D'^ > involving only spinless particles is dimensionless. 
In the {E\,E2) plane, the phase space is defined by the constraints 

< Ex < 



with 



2 niDs 

E_{Ei) < E2 < E+{Ei) , (4) 



E±{E) = -{mD, - E) ± - ^ 2 ^ 2 o 77- \ ■ 5 

2' " '21 mjj +m^-2mD,E J 



Of course, the mass difference between charged and neutral pions being neglected, we have the 
same phase space in the two other planes : {E2,E^) and (EijE^). 

At the quark level, the decay Df Svr is described by a W annihilation diagram, and the 
transition matrix element is given by 

< 37r|r|D+ > = ^ Vud ai < MA^\0 > . (6) 

where oi is the phenomenological parameter introduced by Bauer, Stech, Wirbel and /d^ is the 
leptonic decay constant of the Ds meson. The matrix element of the divergence of the weak axial 
current between the vacuum and three pion final state is an unknown form factor F(Ei, E2), function 
of two independent variables, taken as the pion energies. 

<37r|A^|0> = niD, F{Ei,E2) . (7) 

Since the transition matrix element is dimensionless, the form factor F{Ei,E2) i is also dimensionless. 
In the (-El, E2) plane, the Dalitz plot is given by : 

dr{Dt ^37r) = ^ (^^J \V,s\' iKdP al fl^ im,i?2)P dE^ dE2 . (8) 

The decay rate is obtained by integration of the distribution (P) over the energies Ei and E2 inside 
the phase space described in Eq.(^). The branching ratio can be written in the form : 



Br{D+ 37r) 



n 64 7r3 



Vcs? \Vu,? al I . (9) 



V2 

where the dimensionless integral / is defined by : 

I=-^[ [ \F{EuE2)\^ dEi dE2 . (10) 

Using = 4.67 • IQ-^^^ Q and /d, = 280MeV i, we get 

Br{Dt Svr) = 1.31 • 10"^ aj I . (11) 

From the analysis of colour favoured D meson decay, a reasonable value for oi is oi = 1.26 |^ and 
we obtain 

Br{Df 3tt) = 2.08 • 10"^ / . (12) 



tQur F{Ei,E2) is related to the F4 of the Ref.|[ by F{Ei,E2) = mo, Fi{si,S2,Q'^ = m|,^ 
^The leptonic decay fi'^ + has been measured by two groups and the extracted values of fo^ are fo^ = 

(232 ±69) MeV and fo, = (344±76) MeV. Our choice of fn, = 280 MeV is consistent with data and with theoretical 

expectations. 
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We make the following choice of vr-meson energy variables ; 

i) final state vr+vr+vr- Ei{7r+), £2(71+), £3(71-) , 

ii) final state 7r°7r07r+ Ei{7r°), E2{tt^), £3(71+) . 

Due to Bose-Einstein symmetry, the function F{Ei,E2) is symmetrical in the exchange between 
El and E2- Such a property extends to the Dalitz plot given by the double differential distribution 
represented by \F{Ei, E2)\'^ ■ It is probably premature to discuss the details of the Dalitz plot and it 
might be interesting to define one meson energy distributions : 

G(Ei) = / \F{EuE2)\^ dE2 , (13) 

rriDs JE.{Ei) 

H{Es) = J_ f''^^''''^ \F{mD,-E2-Es,E2)\^ dE2 , (14) 
rriDs Je.(Ez) 

where G{Ei) corresponds to the 7r+(7r'^) for the final state 7r+7r+7r^(7r''7r'^7r+), and H[E-i) corresponds 
to the 7r^(7r+) for the final state 7r+7r+7r^(7r'^7r'^7r+). 

The presently available data concern the rates for the decay mode Df — > vr+Tr+vr^. They are given 
in Table 1, the last column indicates the value of the quantity / determined from experiment |||] by 
using Eq.(|l2D. 



Mode 


Experimental Branching Ratios 


Experimental values of / 




(1.35 ±0.31) • 10-2 


0.649 ±0.149 


(7r+7r+7r")Ar/j 


(1.01 ±0.35) • 10-2 


0.486 ±0.168 


/ovr+ 


(1±0.4) • 10-2 


0.481 ±0.192 




< 0.28 • 10-2 


< 0.135 



Table 1 



III. Phenomenological Model for Stt 

In order to describe the decay mode Df 3tt, we use a simple phemonenological model where 
the function F{Ei, E2) is the sum of two contributions ; 

F{Ei,E2) = Fnr + Fres{Ei,E2) . (15) 

Obviously, this decomposition is guided by experimental data . The first component F/vj? is assumed 
to be a real constant describing the non resonant part of the amplitude. The second component 
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Fres{Ei, E2) is associated to the quasi two body final state — > /o + tt^ followed by the subsequent 
decays /o — > vr+vr" or vt'^tt'^. 

1). At first, consider the charged case /o — > vr+vr^. The vr" meson has an energy and the two 
vr"'' mesons energies Ei and £'2- We then have two possible resonant combinations and the general 
form of Fres{Ei,E2) is : 

Fres{EuE2) = D, {BW{Ei) + BW{E2)} , (16) 

where the dimensionless Breit-Wigner function BW{Ej) is defined by : 

2 

rrif 

BW{E,) = ^° ^ (17) 



m 



fo 



with Sj and Ej related by Sj = m|)^ + — Itud^Ej { j = 1, 2). 

The complex constant Dc corresponds to the transition W'^ foT^^, for which we assume, by the 
Partial Conservation of the Axial Current {PC AC) in Eqs.(^) and (|^), the existence of an intermediate 
state having the quantum number of a vr"^ meson ^, ^ and which might be the vr meson itself or its 
recurrence 7r(1300). It turns out that the vr meson intermediate state will give a contribution many 
order of magnitude smaller than that of the 7r(1300) and only the later one is retained given the 
following expression for Dc : 

-^c = ^ — f T — g^' f 9foTT+iT- ■ (18) 

In Eq.(|T^), vr' = 7r(1300), and the dimensionless coupling constants 5M0M1M2 describe the decay 
of a spinless meson Mq into two spinless mesons Mi and M2 ■ 

< Ml M2 I T I Mo > = mo 9MoMiM2 ■ (19) 

The term m^i^ in the numerator of Eq.(|l8|) comes from the definition Eq.(|^), the term in the 
denominator of Eq.(^) when combined with Eq.([T7|) yields the definition Eq.(|l9|). Also m^i comes 
from Eq.(l^) and f^i comes from the coupling between vr and W. The numerical value of (7M0M1M2 
is obtained from the experimental decay rate r(Mo Mi + M2) by the relation 

2 r(Mo^Mi + M2) 
5moMiM2 = — , (20) 

where K is the final momentum in the Mq rest frame. 



^We have used for width the simple energy dependence r(s) = -^r(m). In the case of the 7r(1300) width, more 
sophisticated dependence has been proposed Because of the not too large difference between the 7r(1300) and Df 
masses, we expect the sensitivity of Eq.(p^) to different forms of F^i (s) to be relatively modest. 
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Now we consider first the case /o — > 27r. Using the experimental parameters 



= 980 MeV , Tf^ = (49 ± 9) MeV , 
Br{fo 27r) = 0.781 ± 0.024 , 



(21) 



we obtain 

= 1.144 ±0.109 , 5j(,^(vo = 0.809 ± 0.077 . (22) 

The isoscalar property of the /o meson has been taken into account for relating both couphng constants. 
Secondly, for the decay 7r(1300) — > /o + tt, experimental data are poor ^ and we shall use : 



The result is 



m / = 1300 MeV 

TV 

Br{n' foTT) = 0.68 

9 -k' fair 



r^/ = (400 ± 200) MeV , 



5.2 



1.2 
1.5 



(23) 



(24) 



The last parameter entering in Eq. (^) is the leptonic constant f^' . The estimate of reference is i 

4, = 40 MeV (25) 

Now we are in a position to compute the numerical value of the dimensionless coupling constant 
Dc written in the form Dc = \Dc\ exp{i ^d)- Retaining only the large error due to the poor knowledge 
of the total 7r(1300) width, we obtain for different F^/ values : 



(26) 



2). In the second case /o — > vr'^vr'', we call Ei and E2 the energies of the vr'^'s and the energy 
of the vr"*". The general form of FjiEs{Ei, E2) is 

Fres{Ei,E2) = Dn BW{E^) , (27) 

and because of Eq. (|2^) , = Dc/V^. 

IV. Results of the model for Tr+Tr+vr" 

The results of the model for the one pion energy distributions and for the rates are now given in 
the three following cases : 





600 




0.3068 




152° 




400 


MeV , \Dc\ = 


0.2679 


, $D = 


160° 




200 




0.1981 




170° 



''Different values of f^i have been proposed in the literature. See for instance references |g| and 
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1) non resonant tt^tt^tt^ state : Fatr 

2) quasi two body /ovr"'' 7r+7r^7r+ state : -Fres 

3) superposition of the non resonant and resonant amplitudes : F/v/j + Fres 

1). In the non resonant case, we have an uniform Dahtz plot and for a constant -F/v_r, we get from 
Eqs.dl) and : 

G^niE) - H^niE) - ^ | ml^+ml-2mn.E 1 " ^'^^ 

Here E stands for Ei, E2, E^ indifferently. 

The energy distribution is represented on Figure 1 for Fj^r = 1. The integral Inr given by 

2 2 

1 



Inr = / ' Gnr{E) dE , (29) 

has the numerical value 

Inr = 0.1053 F^R . (30) 
Using, for Inr, the experimental value quoted in Table 1 in the non resonant case, we obtain 

IF^^I = 2.15 t Ut ■ (31) 

2). Now we consider the quasi two body case fo + t^~^ Tr^TT~TT~^ for which the amplitude 

Fres is given by Eq.(16). It is convenient to write 



\FREsiEi,E2)\^ = \D,\^ K,{Ei,E2) , (32) 
where the function Kc{Ei, E2) is the sum of three terms ; 

Kc{Ei,E2) = \BW{Ei)\^ + \BW{E2)\^ + 2Re{BW{Ei)BW{E2y} . (33) 



After integration of Kc{Ei, E2) over E2 at fixed Ei, as indicated in Eq.(13), we obtain the vr"*" meson 
energy distribution in the form : 

Gres{Ei) = \D,\^ Kj,+\Ei) . (34) 

The function Kc^\Ei) is represented on Figure 2. Clearly we see the narrow peak due to the Breit- 
Wigner term in Ei and a quasi-constant background corresponding to the Breit-Wigner term in £"2. 
The interference between the two Breit-Wigner terms gives a very small contribution to Kc'^\Ei). 



**In the limit rrin ~ 0, the calculation of Inr is trivial and the result is Inr ~ 0.125 Correction due to nin 7^ 

is large and of the order of 16% in spite of the very small value of m^/mn^ ~ 0.07 . 
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Integrating now Kc{Ei, E2) over E2 at fixed E^, as indicated in Eq.(|l^), we obtain the vr meson 
energy distribution : 

Hres{Ez) = \Dc\^ K^'HEs) . (35) 

The function \e?j) is also represented on Figure 2. The result is a quasi-constant plateau in most 
of the allowed phase space domain of E^,. 

In order to obtain the rate, we must perform a second energy integration. The quantity Ires is 
written in the form : 

Ires = \Dc\^ , (36) 

where Kc is given by : 



K, = ^ K^-^KE,) dE, = ^'""^ K^iE,) dEs . (37) 



t\Ei) dEi = ^ I ^<^(-) 
The numerical value obtained for Kr is Kr = 5.7217. 



Using now the results of Eq.(pq) for |-Dc|i we obtain the numerical value in our model of Ires 
defined in Eq.(^). The result is 

Ires = 0.411 t '^iH , (38) 

where, as previously explained, the errors correspond to F^/ = 400 it 200MeV. 

The theoretical value of Ires in Eq.(38) is consistent with the measured value of 0.481 it 0.192 



quoted in Table 1. If we consider seriously the experimental one standard deviation limit. Ires ^ 
0.289, we get a lower bound for F^/, F^/ > 263MeV (assuming f^/ = 40Mey). There is no upper 
bound constraint on F / . 



3). Finally we consider the full amplitude written in Eq.(15). As the constant Fj^r is real, we get 



\F{Ei,E2)\^ = F^ji + 2 FNRRe{FREs{Ei,E2)} + \Fres{Ei,E2)\^ . (39) 



We integrate over Ei and E2 in order to obtain the quantity / of Eq.(lO). 
Using the previous results Eqs.(|30|) and (|3^), we obtain 



7 = 0.1053 F^R + 2 Fnr \D^\ \J\ C os{^ d + ^j) + ^-7217 \Df , (40) 



where the complex integral J is defined by : 



J=-^ I I {BW{Ei) + BW{E2)} dEi dE2 = \J\ e^*^ . (41) 



The results are 

|J|= 0.2677 , $j = 89.60° . (42) 



Using the experimental constraints for Inh and / given in Table 1, from Eqs.(|30|) and (^0|), we 
can check the consistency of our model by determining the parameter -Fatr and compare to Eq.(ll). 
The results are presented on Table 2 for three value of F^/ , T^/ = 600, 400, 263 MeV. For large 
value of F^/, we obtain only positive solutions for F/v_r. When F^/ decreases, \Dc\ decreases and the 
interference term in Eq.(|40|) becomes smaller, then we obtain both positive and negative solutions for 



F^/ (MeV) 


Fnr 


Inr 


/ 


Ires 


600 


1.7378 


1.9898 


0.318 


0.417 


0.719 


0.798 


0.5385 


400 


1.7378 


2.1680 


0.318 


0.495 


0.642 


0.798 


0.4105 


263 


1.7378 


2.3385 


0.318 


0.576 


0.557 


0.798 


0.2890 




-1.7378 


-2.0671 


0.318 


0.456 


0.657 


0.798 


0.2890 


Experiment 




0.318 - 0.654 


0.500 


- 0.798 


0.289 - 0.673 



Table 2 



The full pion energy distributions coming from the total amplitude in Eq.(p9|) depend on the two 
quantities F/vj? and Dc- Writing the squared modulus of the total amplitude in the form : 

\F{E^,E2)f = F^ji + 2FNR Re{D,MEuE2)} + \D,f Kc{E^,E2) , (43) 

where Kc{Ei, E2) is given by Eq.(^) and Jc{Ei,E2) is simply given by : 

JciEi,E2) = BW{Ei)+BW{E2) . (44) 

Of course, K,{EuE2) = 1/^(^1,^2)!'. 

After integration of Eq.(^3|) over E2 at fixed Ei, we obtain the vr"'' meson energy distribution 
GciEi) : 

Gc (El ) = F^ji Gnr {El) + 2 Fnr Re{Dc J^{Ei)} + \Dc\^ K+{Ei) . (45) 

The quantities Gnr{Ei) and K^{Ei) have been represented on Figures 1 and 2, respectively, and 
J+(£'i) is defined by : 

1 rE+{Ei) 

J+{Ei) = / J,{Ei,E2) dE2 . (46) 

In a similar way, integrating over E2 at fixed E's, we obtain the 7r~ meson energy distribution in 
the form : 

^,(^3) = F^j, Gnr{E3) + 2 Fnr Re {D, J-{E^)} + K-{E^) . (47) 
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In Eq.(|47|), GNRiE^) as given by Eq. (|28|) is represented on Figure 1. The quantity (E^) is repre- 
sented on Figure 2 and J~{E3) is defined by : 

1 rE+{E3) 

J-iEs) = / Je(mD^ -E2- Es, E2) dE2 . (48 

As mentioned above, the full vr^ and vr^ energy distributions Gc{Ei) and Hc{E^) are functions of F^r 
and Dc- They are represented respectively in Figures 3 and 4. 

V. Prediction of the model for -k^tx^-k^ 

We have no experimental data on the decay mode D'^ However our model can make 

predictions for the three following cases : 

1) non resonant vr^vr'^vr"'" state, 

2) quasi two body /o'/r+ — > 7r''7r''7r+ state, 

3) superposition of the non resonant and resonant amplitude. 

1). In the non resonant case of a constant function F]\[ji{Ei, E2), we obtain the same shape for 
one pion energy distribution as shown on Figure 1. 

However the constant Fnr has no reason to be the same for the two modes Tr^vr+vr^ and 
We now discuss this problem. By assumption, with a constant function Fnr{Ei,E2), we have a full 
symmetry in space between the three pions. From the Bose-Einstein symmetry, the isospin config- 
uration has also to be totally symmetric. Consider now a third rank fully symmetric tensor in a 3 
dimensional space. It has 10 independent components. With respect to the isospin S0(3) orthogonal 
group, such a tensor is reducible into an isospin / = 3 part with 7 components and an isospin 7=1 
part with 3 components. In our specific problem, only the later part contributes, the ud weak current 
being an isovector. By inspection of the relevant Clebsch-Gordan coefficients, we obtain the result : 

Fnr{^+^^^-) = 3 FjvK(7rW+) . (49) 

As a consequence, in our model, we obtain for the non resonant part : 



11.1 % , (50) 



T{D+ ^ 7r+7r+7r-)jvi? 

and the non resonant branching ratio Br{Df TT^-K^7r^)NR is expected to occur only at the 10~^ 
level. 

2). For the quasi two body decay of Df foT^^ tt^tt^tt^, the amplitude Fjies is given by 
Eq.dl) : 

\Fres{EuE2)\^ = \Dn\^ Kn{Ei,E2) (51) 
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where the function Kn(Ei, E2) is simply given by : 

Kn{Ei,E2) = \BW{E^)\^ = \BW{mD^-Ei-E2)\^ (52) 

After integration over E2 at fixed Ei, as indicated in Eg. (|l^) , we obtain the vr'^ meson energy distri- 
bution in the form : 

Gres{Ei) = \Dn\'' kP{E,) (53) 

The function K^\Ei) is a quasi flat distribution represented on Figure 5. 

Integrating now over E2 at fixed E^, as indicated in Eq.(^4|), we have the vr"'" energy distribution 
in the form : 

HREsiEs) = \DM^ K^n\Es) (54) 

The function K^\e^) is represented on Figure 5 and the result exhibits a narrow Breit-Wigner peak 
similar to the one drawn on Figure 2 for K^{Ei). 

The decay rate is now obtained by performing a second energy integration. The quantity Ires is 
written in the form : 

Ires = \Dn\^ Kn , (55) 

where K]\j is given by : 

1 



Kn = ^ I K^;^\e,) dE, = / ''""^ kI^\e,) dE, (56) 

The numerical value of Kn is Kn = 2.7469. 

Using the relation between and Dc, we make the prediction 

The departure of this ratio from 25 % is simply due to the interference between the two Breit-Wigner 
contributions in Df — > foir^ vr+vr^Tr^. 

3). We finally consider the full amplitude written in the form of Eq.(|l5|). An equation similar 
to Eq.(]40|) gives the total rate, and going from the mode 7r'*"7r+7r^ to the mode vr^^Tr'^Tr"'' implies the 
changes : 

Fnr^^Fnr , D.^^D, , J^Ij , K.^Km (58) 
and we obtain instead of Eq: (^) : 

40^0^+ = 0.0117F^^ + 0.2357FjvK|I?c||J|Cos(^>D + $j) + 1.3735|Z?eP • (59) 
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Using the solutions for Fjyji given in the Table 2, we can compute the decay rate for the 
mode as a function of -Fatr and of the vr width T^/ via the quantity Dc- In Table 3, we show the 
results for the ratio R defined by : 

r(7r+7r+7r-) ^ ^ 



r^' (MeV) 


Fnr 


R 


600 


1.7378 1.9898 


20.2 ±0.5 % 


400 


1.7378 2.1680 


18.5 ±0.8 % 


263 


1.7378 2.3385 
-1.7378 -2.0671 


16.7 ±1 % 
16.3 ±0.5 % 



Table 3 



The pion energy distributions depend on Fj\fR and Dc, and we have for the squared amplitude an 
expression similar to Eq.(^) with the appropriated changes indicated on Eq.(^) : 



l^NR + ^FNRRe{DcJN{EuE2)} + ^\Dc\^Kn{Ei,E2) , (61) 



where Jn{Ei, E2) contains only one Breit-Wigner amplitude : 

Jn{Ei,E2) = BWiEs) = BW{mD,- El- E2) , (62) 

and i^TV (£^1,^2) = |J7v(^i,^2)|'. 

The full vr'^ energy distribution is obtained after integration of \F{Ei, E2)\'^ over E2 at fixed Ei, 
and the full vr^ energy distribution after intergration over E2 at fixed E^ : 

Gn{Ei) = \Flj,GNR{Ei) + ^FNRRe[Dcjf{Ei)} + \\Dc\^K^^\Ei) , (63) 
Hn{E3) = If^nrGnr{Es) + ^F^RRe{Djj^\Es)} + ^\Dc\^kI^\E3) , (64) 

The functions K^\Ei) and K'^{E^) have been represented on Figure 5. The functions and 
are defined by : 



jPiEi] 



J^n\Es) 



1 rE+iEi) 



/ BWiniD, -El- E2) dE2 , 

Je^{Ei) 

BW{E^) dE2 , 



1 fE+iEa) 



(65) 
(66) 
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The TT^ and vr^ energy distributions Gn{Ei) and Hn{E^) are functions of F^r and Dc- They are 
represented on Figures 6 and 7 with the same values for F^/ and -F/v_r as those used in Figures 3 and 4. 
The first observation is the minor role now played by the constant non resonant component essentially 
due to the isospin factor Eq.(^9|). However the sign of the interference between the two components 
can be observed at large values of as illustrated in Figures 7-a and 7-b. 

V. Summary and Concluding Remarks 

The large value experimentally observed for the branching ratio of the decay mode D'^ — > vr+vr+Tr" 
is a very interesting problem which necessitates a better understanding of the role played by the 
W annihilation mechanism in D meson decay, and which is probably more important than usually 
expected. 

We have proposed to construct the unique structure function F{Ei,E2) as the sum of a constant 
non resonant term described by a real parameter -F/v_r, and a Breit-Wigner resonant term associated 
to a quasi two body /o(980) + vr"^ state which seems to play an important role as indicated by 
experiments. The quasi two body + tt state which is commonly considered by previous authors 
IQ, ^ |6| is experimentally depressed and has been disregarded here. 

It is possible to determine the resonant amplitude using parameters experimentally known, even if 
a large uncertainty remains on their precise experimental values. We have only retained, for simplicity, 
the uncertainty due to the 7r(1300) width. Our model depends only on one free parameter, and we 
show that it is possible to find values of -F/vj? fitting the three observed rates : total Df vr^vr+Tr^, 
non resonant {Df — > 7r+7r+7r^)Arij, and resonant Df — > /ovr+ Tr+vr+vr^. This result is obviously 
non trivial. 

Of course, the Dalitz plots are determined by the squared modulus of the structure function 
\F(Ei, E2)\'^ and they are predicted by our model. Before the double differential quantities could be 
observed for the reason of statistics, it would be easier to measure the one pion energy distributions 
which are nothing but projections of the Dalitz plots on one energy axis. The figures shown in this 
paper present the various situations to be compared with experimental results when available. 

The decay mode Df — > 7r'^7r^7r+ has not been experimentally observed, however the rates, the pion 
energy distributions, and the Dalitz plots can be computed in our model. We compare the vr^vr'^vr^ 
and Tr+vr+vr^ integrated widths through the ratio R in Eq.(|60|) and Table 3, and we also draw the vr'^ 
and 7r+ meson energy distributions. The quantity R is predicted to be in the 16% — 20% range and 
that might be the reason why the mode has not yet been detected, its branching ratio 

being expected to be few 10~^. 
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We wish to emphasize that in the absence of additional experimental informations, our crude model 
is the simplest one accounting for the observed rates. In particular, for the resonant component, the 
normalization is estimated theoretically and constrained experimentally. Fortunately both approaches 
overlap, a fact a priori not evident, and they can be improved when experimental uncertainty will be 
reduced. 
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Figure captions 



1. Figure 1 : The shape of the pion energy distributions for a uniform Dahtz plot. E stands 
for El, E2, E3 indifferently. 

2. Figure 2 : The shapes of the 7r+ and tt~ meson energy distributions for Df — > /o7r+ — 
7r+7r~7r+. K^{E) represents the energy distribution of 7r+ meson and K~{E) represents that of 
7r~ meson. The quantity Eq = 7.45 GeV is associated to the /o resonance. 

3. Figure 3 : The fully normalized vr"'" meson energy distribution for D'^ — > 7r"'"7r~7r"'". 

3-a) for r^/ = 600 MeF, 1.74 < Fnr < 1.99 

3-b) for r^/ = 263 MeV, -2.07 < Fnr < -1.74 

The two lines in each curve correspond to the extremum values of Fnr. The quantity £^0 = 
7.45 GeV is associated to the /o resonance. 

4. Figure 4 : The fully normalized tt^ meson energy distribution for Df — >■ 7r'''7r~7r'''. 
The parameters are the same as in Figure 3. 

5. Figure 5 : The shapes of the 7r° and tt"^ meson energy distributions for /o^^ 

The quantity Eq = 7.45 GeV is associated to the /o resonance. 

6. Figure 6 : The fully normalized 7r° meson energy distribution for Df — 7r°7r°7r+. 
The parameters are the same as in Figure 3. 

7. Figure 7 : The fully normalized 7r+ meson energy distribution for Df —>■ Tr^Tr'^Tr"'". 

The parameters are the same as in Figure 3. The quantity Eq = 7.45 GeV is associated to the 
/o resonance. 
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Figure 1: The shape of the pion energy distributions for a uniform Dahtz plot. E stands for £^2, 
indifferently. 
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Figure 2: The shapes of the 7r+ and tt" meson energy distributions for Df /o7r+ 7r"'"7r~7r"'". 
K^{E) represents the energy distributions of tt"*" meson, and K^{E) represents that of tt^ meson. 
The quantity Eq = 7.45 GeV is associated to the /q resonance. 
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Figure 3: The fully normalized 7r+ meson energy distributions for — >■ 7r+7r^7r+. 3-a) for V^i = 
600 MeF, 1.74 < Fnr < 1.99, 3-b) for T^- = 263 MeV, -2.07 < Fnr < -1.74. The two lines in 
each curve correspond to the extremum values of F^r- The quantity Eq = 7.45 GeV is associated to 
the /o resonance. 
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Figure 4: The fully normalized tt^ meson energy distributions for D+ Tr+Tr^Tr"*". The parameters 
are the same as in Figure 3. 
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Figure 5: The shapes of the 7r° and tt"*" meson energy distributions for D+ /oTt"*" 7r°7r°7r+. The 
quantity Eq = 7.45 GeV is associated to the /□ resonance. 
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Figure 6: The fully normalized 7r° meson energy distributions for Z)+ Tr^Tr^Tr"*". The parameters 
are the same as in Figure 3. 
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Figure 7: The fully normalized tt"*" meson energy distribution for Z)+ Tr^Tr^Tr"*". The parameters are 
the same as in Figure 3. The quantity Eq = 7.45 GeV is associated to the /o resonance. 
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